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Recent paper by Nataraj et al (Phys. Rev. Lett. 106, 200403 (2011)) presents calculations 
of the EDM enhancement factor for Tl, which disagrees with previous most accurate calculations. 
The authors claim that their calculations of Tl EDM are the most accurate due to more complete 
treatment of higher-order correlations. In this note we argue that this claim is not supported by 
sufficient evidence. Nataraj et al also present misleading comments about our calculations. We 
explain our method and reply to the Nataraj et al comments. 
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Recent paper [l| presents calculations of the EDM en- 
hancement factor for Tl, which disagrees with previous 
most accurate calculations 0, The calculations of 
this kind are used for a search for new physics beyond 
the standard model in atomic experiments. Therefore, 
it is important to use the most accurate results for the 
interpretation of the experiments. 

The authors of [l| make many misleading comments 
about our calculations of thallium EDM [3j. Contrary 
to the statement in [l[ the atomic electric field interact- 
ing with electron EDM is calculated in Q as a deriva- 
tive of the total potential which includes both nuclear 
and electron parts. This is done in a same way as in 
our early work on EDM of Fr and Au [3] ■ The formula 
E = Zer/r 3 for the leading contribution to the atomic 
electric field presented on first page of [3| may indeed 
make an impression that only nuclear field is included. 
However, the formula (2) few lines below clearly includes 
screening functions Q(r) and P(r) for both the nuclear 
Coulomb and the external electric field. By the way, the 
inclusion of the electron electric field change the matrix 
elements of the electron EDM for thallium by 0.4% only. 
This is because the main contribution comes from short 
distances where the electron electric field is small since 
the electron potential rapidly tends to a constant inside 
the Is orbital. 

The authors of [l| claim that atomic core is strongly 
contracted in the V N ~ 3 zero approximation used in our 
calculations [3]. In fact, it is not. Fig [1] shows the out- 
ermost 5^5/2 core function of Tl calculated in V N ~ 3 and 
V^" 1 approximations. The difference between the func- 
tions is very small. This is due to the fact that the va- 
lence 6s and 6p electrons are located outside of atomic 
core. Their charge distribution creates almost constant 
potential and no electric field inside the sphere where all 
inner electrons are located. Therefore, the valence elec- 
trons have practically no effect on the core wave func- 
tions (see [5[ for a detailed discussion). The change is 
even smaller for other core functions. The core functions 
enter the configuration interaction (CI) Hamiltonian via 
core potential V core to which all core electrons contribute 
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FIG. 1: The 5d 5/2 core function of Tl calculated in F" 3 
(solid line) and V N_1 (dashed line) approximations. 

(from Is to 5d). The difference for V core in V N ~ 3 and 
V N ~ X approximations is very small Moreover, the 
corresponding corrections to the configuration interac- 
tion (CI) Hamiltonian have been included in [3J using 
the many-body perturbation theory approach. 

As it is well known, the eigenstates of a Hamiltonian do 
not depend on the basis one uses. The valence states are 
indeed different in the V N ~ 3 and V N ~ 1 approximations. 
However, this should have no effect on the final results 
as long as the saturation of the basis for valence states is 
achieved. There are only two conditions the basis states 
must satisfy: (a) they must be orthogonal to the core, 
and (b) they must constitute a complete set of states. 
Therefore, it does not matter whether valence states are 
calculated in the V N ~ 3 or V*" 1 potential or by any other 
method (e.g., a popular basis ip n (r) = r n ip (r) @), the 
final results should be the same. If there is any differ- 
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ence in the results, the most likely reason for this is the 
incompleteness of the basis set. 

In spite of no difference in final results there is a good 
reason for the use of the V N ~ 3 approximation - the sim- 
plicity and good convergence of the many-body pertur- 
bation theory (MBPT) for the core- valence correlations. 
When an approximation different from the V N ~ 3 is used 
one has to include the so called subtraction diagrams [7|, 
while there are no such diagrams in the V N ~ 3 approx- 
imation. Large energy denominators suppress the value 
of the correlation terms in the V^ -3 approximation en- 
suring good convergence of the MBPT [5[. There must 
be large cancelation between subtraction and other dia- 
grams to ensure the same final results if any other ini- 
tial approximation is used. This is very similar to the 
well known fact that the Hartree-Fock basis is the best 
choice for any MBPT calculations. Initial approximation 
might be better in some other approximation, however, 
strong cancellations between the subtraction and other 
diagrams would lead to poor convergence of the MBPT. 

The authors of [l[ claim that the major drawback of 
our work Q is the inclusion of the core- valence correla- 
tions in the second order only. However, the correlations 
between the valence electrons and core electrons below 
the 6s state are small which is evident from the fact that 
their inclusion change the EDM of Tl by 3% only 
Therefore, only the correlations between three valence 
electrons should be treated to all orders. This is done 
in Q to a very high precision using the CI technique. 

Early calculations of thallium EDM by Liu and 
Kelly |2J were performed by the same relativistic coupled- 
cluster (RCC) method as those used in Table I of jl| 



presents term by term comparison between the contribu- 
tions to the enhancement factor calculated in both works. 
For some terms the agreement is perfect, for others there 
is strong disagreement. The authors of 1] claim that 
the disagreement is due to more accurate treatment of 
higher-order correlations in their work. It might be pos- 
sible to prove this claim by switching off the higher-order 
terms and reproducing the results of |2|. The benefit of 
having the test is enormous. Without it, no other reasons 
for disagreement can be excluded (e.g., incompleteness of 
the basis) . There is no indication that such test has been 
performed in 

In conclusion, we would like to note that the calcula- 
tions of the Tl EDM due to electron EDM [l[ and SPS 
interaction [8J do not satisfy a simple consistency test: 
the ratio of the EDMs due to two operators must be ap- 
proximately equal to the ratio of the s — p single-electron 
matrix elements of these operators. This ratio is the same 
for all important single-electron matrix elements. This is 
because only short distances, where single-electron ener- 
gies can be neglected, contribute to the single-electron 
matrix elements of the T,P-odd operators. The ratio ap- 
proximately equals to 89d e /C SP 10~ 18 e cm according to 
analytical estimates, or 83d e /C SP 10 -18 e cm in more ac- 
curate numerical calculations [9j. This ratio must hold 
in any order of the MBPT if the s 1/2 - Px/2 T,P-odd 
matrix elements dominate at the Hartree-Fock level (see 
for details) . The ratio of the results of [l[ and Q is 
115d e /C 10 -18 e cm. This may indicate that impor- 
tant many-body effects are missed in one (or both) of the 
works 0,11. 
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